Localized mode interactions in O-ir Josephson junctions 



Hadi Susanto 1 and Gianne Derks 2 

'School of Mathematical Sciences, University of Nottingham, University Park, Nottingham, NG7 2RD, UK 
2 Department of Mathematics, University of Surrey, Guildford, Surrey, GU2 7XH, UK 

A long Josephson junction containing regions with a phase shift of tt is considered. By exploiting 
the defect modes due to the discontinuities present in the system, it is shown that Josephson 
junctions with phase-shift can be an ideal setting for studying localized mode interactions. A 
phase-shift configuration acting as a double-well potential is considered and shown to admit mode 
tunnelings between the wells. When the phase-shift configuration is periodic, it is shown that 
localized excitations forming bright and dark solitons can be created. Multi-mode approximations 
are derived confirming the numerical results. 
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Introduction. A Josephson junction is a system con- 
sisting of two layers of superconductors separated by a 
nonsuperconducting barrier. Electrons forming so-called 
Cooper pairs can tunnel across the resistive barrier even 
when there is no applied voltage difference. Theoretically 
predicted by Josephson^ and first observed experimen- 
tally in Ref. 0, the only requirement for the occurrence of 
Josephson tunneling is a weak coupling of the wave func- 
tions of the two superconductors. The supercurrent (I s ) 
is proportional to the sine of the electron phase-difference 
across the insulator (u), i.e. I s ~ sin^. 

Bulaevskii et al. 3,4 proposed that a shift of tt can be 
introduced in the phase difference u of a Josephson junc- 
tion by installing magnetic impurities, which has been 
confirmed recently 5 -. Present technological advances can 
also impose a 7r-phase-shift in a long Josephson junction 
using various means, including multilayer junctions with 
controlled thicknesses over the insulating barrier 6,7 , pairs 
of current injectors^, and junctions with unconventional 
order parameter symmetry^"—. 

New phenomena may occur when a junction with 
phase-shifts is connected to a normal junction, i.e. O-tt 
Josephson junctions. These include the presence of a 
half magnetic flux quantum induced by spontaneously 
created supercurrent circulating in a loop. Such unique 
characteristics offer promising future device applications, 
such as novel circuits for information storage and pro- 
cessing in both classical and quantum limits 12 and arti- 
ficial crystals for simulating and studying energy levels 
and band structures in large systems of spins 13 (see also 
Ref. and references therein). Here, we demonstrate 
that Josephson junctions with phase-shifts is an ideal set- 
ting for showcasing many interesting features of localized 
mode interactions. Arguably the dynamics of the super- 
conductor phase-difference can be seen to be analogues to 
that of atomic wave functions of Bose-Einstein conden- 
sates (BEC)^ 5 -"— in an external potential. In particular, 
we consider Josephson junctions with phase-shift config- 
urations acting as a double well and a periodic potential. 

The interesting phenomenon of mode tunneling in 
BECs in a double well potential was predicted by Smerzi 
et al. 20,21 , followed by experimental observations 2 ^ 2 ^. 
Here, we will show that defect modes due to the pres- 



ence of phase-discontinuities in O-tt Josephson junctions 
can be exploited to observe a similar mode tunneling, 
which can also be viewed as Rabi oscillations of two in- 
teracting modest Periodic defects exhibiting mode self- 
trapping analogues to BECs in optical lattices 25 will also 
be discussed. 

Double-well potential. A O-tt Josephson junction with 
the superconductor phase difference u at position x and 
time t is described by the sine- Gordon equation 



= —6{x) sin?/, 



(i) 



where x and t have been normalized to the Josephson 
penetration depth Aj and the inverse plasma frequency 
respectively. The function 0{x) is piecewise con- 



stant representing the presence of tt junctions. A double 
well potential with two tt-j unctions of length a separated 
by a 0-j unction with length 2L is described by 



0{x) 



{-:; 



L < \x\ < L - 
elsewhere. 



(2) 



At the points of discontinuities, the boundary conditions 
are 

lim u(x.t) = lim u(x,t): 

x^>±{L,L+a} + x^±{L,L+a}~ ,„\ 

lim u x (x,i) = lim u x (x,i). ^ ' 

x^±{L,L+a} + x^±{L,L+a}~ 

If (j) solves (pQ), the linear stability of the solution can 
be analyzed by substituting the spectral ansatz u = (j) + 
v(x)e xt and linearizing about |H|oo small to yield the 
eigenvalue problem v xx — X 2 v = 6cos((j))v. 

Equation (pQ) has two constant solutions (mod 27r), 
u = and u = tt. The solution u = tt has unstable 
continuous spectrum and hence is always unstable^ The 
solution u = has stable continuous spectrum A 2 < — 1 
and the discrete spectrum (eigenvalues) can be calculated 
analytically. 26 Indeed, the largest two eigenvalues A± of 
u = solve the equation 



1 



A 2 ± 



tan ^yjl - A±aJ 



+ A± T e"V 1+A ± l = 0. (4) 
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The corresponding eigenmodes <&±(x) of A± ar< 



,26 



e - y /l+A±(x-L-a) 

cos(y / T^/\±(x — L - 
+ C sin( v /l - A ± (x ■ 



a)) 
- L - 



■ A±x) 



if + cosh(Y / T _ 
+ sinh(^/TTA±x), 



x > L + a; 



a)), L < x < L + a; 



< x < L, 



where C 



i±A± and ± = 



(5) 

2e"\ /1+A ±- L sin(y'l-A ± a) 

As the linearisation operator is a Sturm-Liouville oper- 
ator and even in x, the eigenmodes are simple and the 
eigenfunction <I> + is an even function and <£_ is odd. The 
two eigenvalues A± are depicted as a function of a for 
fixed L = 2 (to the left of the vertical dashed line). 

It is clear that u = has a stability window. The 
change of stability occurs at a critical distance a c when 
the critical eigenvalue crosses the horizontal axis A = 0. 
Using the expression in ((4]), it follows that the critical 
length a c (L) is a c = arctan(e 2L ). 



The translations x<i are determined by the boundary con- 
ditions ([3]) . The non-uniform ground state and its eigen- 
values are presented in Fig. Q] (to the right of the vertical 
dashed line). 

Mode tunneling. In the following, let us first consider 
L = 2a = 2. For those parameter values, u = is a 
stable ground state. The numerically obtained time dy- 
namics of an initially localized excitation in the left well 
is presented in the top panels of Fig. [2j clearly showing 
mode tunneling. Compare it with the time dynamics of 
BECs reported in Refs. HHH 



(a) 





FIG. 1: The eigenvalues of the ground state as a function of a 
for L — 2. The dashed vertical line indicates the bifurcation 
point of the non uniform ground state, where on the left and 
on the right of the vertical line u — is stable and unstable, 
respectively. The dashed lines show the eigenvalues of the 
uniform solution in its instability region. The inset presents 
a ground state when L — 2 and a = 1.65. 

In the instability region, a non uniform time- 
independent sign-definite ground state ±u nc (x) bifur- 
cates. Its expression can be written in terms of Jacobian 
elliptic function s 26 ! 27 as 

{4 arctan(e _a;+:Co ), x > L + a; 
2 arcsin(mi sn(x — xi, mi)), L < x < L + a; 
7r + 2 arcsin(m2 sn(x — #2, ^2)), < x < L. 

(6) 

The parameters mi and m2 are linked to the lengths a 
and L by 

a — 2K(m\) — am _1 (7r/4, mi) 

— am _1 (arcsin(Y // 2(l + m 2 — m 2 )/2mi), mi) 
L = K(rri2) — am _1 (arcsin(Y // 2(l + m 2 , — m 2 )/2m2), m^) 





FIG. 2: (Color online) Panel (a) shows the time dynamics of 
an initially localized excitation in the left well with L = 2 and 
a = 1. In panel (b), the oscillation amplitudes of the dynam- 
ics presented in panel (a) are plotted against time. Panel (c) 
depicts the oscillation amplitude on top of a nonzero back- 
ground with a small initial amplitude for a double well with 
L = 2 and a = 1.65. Coupling coefficients appearing in © as 
functions of a for L — 2 are shown in panel (d). In calculat- 
ing the coefficients, the eigenfunctions have been normalized 
to 1 1 <I>± 1 1 00 = 1. Approximations obtained from the two- mode 
equation J8} are shown as red dashed curves. 

Next, we consider a parameter combination of L = 2 
and a = 1.65, i.e. u = is unstable. In the instability 
region of the constant solution, excitations will oscillate 
on a non-zero background u nc ((6]). The oscillation am- 
plitude in both wells as a function of time is presented 
in panel (c) of Fig. [2l When the initial oscillation am- 
plitude is large enough, we interestingly obtain chaotic 
oscillations (not shown here). 

Two-mode approximations. We will explain the ob- 
served mode tunneling using a two-mode approximation. 
Looking for the solution of the time-dependent equation 
(PQ) of the form 

u(x, t) = u gs + + , (7) 

where u gs is the ground state of the system, i.e. u gs = 



3 



and u gs = u nc when a is respectively on the left and right 
of the vertical line in Fig. [TJ substituting the ansatz ([7]) 
into (pQ), and projecting the equation onto <&± will yield 
up to 0(A n B 4 - n ), n = 0, . . . , 4 



A = 
B = 

with 



A+A + (dA 2 + C 2 B 2 ) + (C 3+ A 3 
A_£+ (C 3 -£ 3 + C 4 -A 2 £); 



C 4+ AB 2 



(8) 



Ci = I 0sm(u g8 )<b\dx, 

^ J -oo 

C 2 = — ± / 0sm(u gs )$+$ 2 _dx, 

^ J -oo 

K± f°° 

^ J — oo 

/oo 
6cos(u gs )§ 2 +$ 2 _ dx, 
-oo 



K± 
6 



(9) 
(10) 

(11) 

(12) 



and K± = ^ &±dx^j . The constants Cj are plotted 

against a in panel (d) of Fig. [2] The internal oscillation 
amplitude u (=t(L + a/2), t) is respectively approximated 
bytx^ + (A(t)±B(t))/2. 

For the uniform ground state (i^ = 0), we have solved 
JHJ) numerically and compared it with the oscillation am- 
plitude of the original equation in panel (b) of Fig.[2j One 
can observe that quantitative agreements are obtained. 
An agreement is also obtained for mode tunneling on a 
non-uniform background, as shown in panel (c) of Fig. 
[2j provided that the tunneling mode amplitude is small 
enough. 

When A(0) and B(0) are large, it is interesting to note 
that even though our two-mode approximation does not 
quantitatively capture the dynamics of the chaotic tun- 
nelings it captures the qualitative transition to chaotic 
behavior. 

Periodic defects. One can also include more phase 
shifts and derive a multi-mode approximation. Below 
we consider the case of periodic shifts alternating n- 
junctions of length a with 0-junctions of length L, i.e., 



1, 



x e k(L + a) + (- 
elsewhere. 



§,§),fcez 



(13) 



It is known that in the limit L — >• oo (i.e., one well), 
u = is stable for a < 7r/4 and unstable otherwise 26 . The 
eigenfunction corresponding to the critical eigenvalue Ao 
of the ground state u gs in the limit L — > oo will be de- 
noted by $>o(x). For a < 7r/4 and u gs = 0, 

-VT+7^(\x\-a/2) 



$oO) 



COs(a/T" 
COs(a/T" 



Aoa)e 



A x) 



< a/2. 



\x\ > a/2, 
(14) 

When L ^> 1, a tight-binding approximation can be 
used to describe the interaction of the defect modes in 
the system with periodic defects. We write 



u(x, t) 



oo 

E 

i=— oo 



A n (t)$ n (x), 



(15) 
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FIG. 3: (Color online) The left panels show the profile of 
numerically exact bright and dark lattice solitons obtained 
from the lattice equation (fT6|) . The insets depict the Floquet 
multipliers of the solution. The right panels present the cor- 
responding time dynamics of (pQ) using the initial condition 
([T5]) with A n (0) shown in the adjacent left panel. All panels 
have L = 10, the upper two have a = 0.5, the lower a — 1.6. 



where <& n (x) — $0 [x — n(a + L)). Performing the same 
procedure as before, one will obtain the lattice equation 

A n = A A n + K 2 A 2 n + K 3 A 3 n 

+ EQKtfij - Ao)A n+J - K 2 A 2 n+J - K 3 Al +j m 
i=±i 

where 

/oo 
<f> (x)$ ±1 (x)dx, 
-00 

/OO 
&o(x) {dl ~ cos u gs ) $ ±1 (:r) dx, 
-OO 

/OO 1 
-sm(u gs )$ 3 (x)dx, 
-OO ^ 

/OO 1 
-C08(Ugs)$$(x) dx, 
-00 ^ 



^2 

K 3 



4 



and Jf(°) = (/^ D $o • Neglecting the nonlinear 

couplings to the neighboring sites {CKi, CK3), the dis- 
crete equation above becomes the lattice equation con- 
sidered by Kivshar 28 , admitting many types of localized 
excitations. In the following, we consider the special 
type, namely unstaggered bright and staggered dark lat- 
tice solitons. In particular, we will show that the lat- 
tice equation ([T6]) can predict the stability of the soli- 
tons in the original equation (pQ) , provided that L is large 
enough. In doing so, we first solve Eq. (fTSj) numerically 
for localized modes using a shooting method and corre- 
spondingly study their stability (see, e.g., the review Ref. 
I29I ) and then use the ansatz (|T5|) at t = as an initial 
condition for the governing equation (pQ). In the follow- 
ing, periodic boundary conditions are used, which are 
relevant experimentally. 30 Examples are shown in Fig. 
[3j Presented in the left and right panels are numerically 
exact solutions obtained from the lattice equation ([T6]) 
and their corresponding time evolution in the original 
equation (pQ). The insets on the left panels depict the 
corresponding Floquet multipliers, where the instability 
is indicated by the presence of eigenvalues lying outside 
the unit circle (dash-dotted line). 

First, we consider the parameter values L = 10 and a — 
0.5, representing the case of stable constant solution u = 
0. Shown in the first and the second row are numerically 
exact bright and dark solitons with the oscillation period 
P = 7.08 (Ao ~ —0.8) and their dynamics. According 
to the lattice equation (fTBjh the bright soliton is stable 
and the dark one unstable. One can note from the right 
panels that the prediction provided by the lattice is in 
agreement with the dynamics in the original system. The 
instability of the dark soliton manifests in the form of the 



destruction of the configuration. 

Next, we consider the parameter values L = 10 and 
a = 1.6, which represent the case of nonuniform ground 
state. The localized mode will then oscillate on a nonzero 
background. Shown in the third and fourth row are nu- 
merically exact bright and dark lattice solitons with the 
oscillation period P = 30 (A « -0.0527). 

According to the lattice equation ([T6]h the bright soli- 
ton has the same stability as the case on stable uniform 
ground state, which is confirmed by the time dynamics 
of the original equation. The stability of the soliton for 
the chosen parameter values is not surprising as the sites 
are rather uncoupled. 

As for the dark soliton, it is interesting to note that 
in the present case it is stable, which is also confirmed 
by the dynamics of the full equation. This implies that 
a nonzero background may act as a stabilizer. Moreover, 
it is also important to note that the modes in different 
lattices have different oscillation fre quen cies. The multi- 
frequency breathers discussed in Ref. I3lll32l may therefore 
be potentially observed in experiments. 

Conclusions. We have considered Josephson junctions 
with phase-shifts of 7r. By exploiting the defect modes 
present due to the phase-discontinuities, the system has 
been shown to be an ideal setting for studying mode in- 
teractions. In particular, we have shown that mode tun- 
neling in a double-well potential can be implemented in 
the system and presented the existence and stability of 
bright and dark solitons in a periodic potential. We have 
shown that the analysis proved by a multi-mode approx- 
imation gives a quantitative agreement with dynamics of 
the original system. 
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